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ABSTRACT 

Using data on the Hi column density distribution in the local Universe, 
fiNni), in this paper we show how to determine g^Nn), the distribution of 
the total gas (Hl+Hll) column density. A simple power law fit to /{Nhi) fails 
due to bendings in the distributions when Nhi < 10^° cm~^ and H is no longer 
fully neutral. If an ultraviolet background is responsible for the gas ionization, 
and g^Nn) oc N^'^, we find the values of a and of the intensity of the background 
radiation which are compatible with the present data. The best fitting values of 
a, however, depend upon the scaling law of the the gas volume densities with 
Nh and cannot be determined unambiguously. We examine in detail two models: 
one in which the average gas volume density decreases steadily with Nh, while in 
the other it stays constant at low column densities. The former model leads to a 
steep power law fit for g^Nff), with a ~ 3.3 ±0.4 and requires an ultraviolet flux 
larger than what the QSOs alone produce at z = 0. For the latter a ~ 1.5 ±0.1 
and a lower ionizing flux is required. The ambiguities about the modelling and 
the resulting steep or shallow A'^^ distribution can be resolved only if new 21-cm 
observations and QSOs Lyman limit absorbers searches will provide more data 
in the Hl-Hll transition region at low redshifts. Using the best fit obtained for 
higher redshift data we outline two possible scenarios for the evolution of gaseous 
structures, compatible with the available data at z ~ 0. 
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1. Introduction 

The Hi column density distribution of intergalactic gas clouds, /{Nhi), is determined 
at high redshifts through the analysis of QSO absorption spectra. The detection of Lyman 
continuum absorption and of Lja absorption lines between us and the QSOs has been used 
to derive /{Nhi) from Nhi ^ lO^^.s ^.^^-2 ^ ;^o2i-5 cm-^ (e.g. Petitjean et al. 1993; 

Rauch 1998; Storrie-Lombardi & Wolfe 2000). Attempts have been made to fit /{Nhi) over 
many decades of Nhi with a simple power law. However, it is more physically meaningful 
to expect a power law behavior for g{NH)- the total (neutral+ionizcd) cohimn density dis- 
tribution. In fact, due to drastic changes in the hydrogen ionization fraction when the gas 
becomes optically thin to photons above 13.6 eV, breaks in /{Nhi) appear for Nhi ~ 10^^- 
10^° cm^^, where a small decrease of Nh corresponds to a rapid change in Nhi- A more 
appropriate approach is to compare the observed /'{Nhi) to the Nhi distribution function 
derived from an assumed g^Nn)-, after ionization changes corrections have been applied. This 
method has been used by Corbelli et al. (2001) to derive the slope of gf(A^jy), assumed to be 
oc Njj"", and the ionization conditions of the gas at redshifts z ~ 2-3. Since the frequency 
of intervening absorption systems is high at these redshifts, it has been possible to limit the 
analysis to high column density absorbers, using data relative to (mostly neutral) Damped 
Lya absorption systems, and to (mostly ionized) Lyman limit systems. A comparison of 
models for dark matter confined systems with the observed /{Nhi) led to values of a > 2 
for z ~ 2-3; an extrapolation of this power law distribution towards systems of lower column 
density, the Lya forest, gives reasonable fits for Nhi > 10^^ cm~^. 

At lower redshifts the number density of absorption systems is much lower and the poor 
statistics based on Lyman limit and Damped Lya alone does not allow a determination of the 
distribution function g{NH)- At zero redshift however these studies can be complemented 
with detailed observations of 21-cm emission line radiation from nearby galaxies. 21-cm 
emission data have the advantage of leaving small uncertainties on Nhi, while saturated 
Lyman limit breaks in QSOs spectra give only lower limits on the intervening Nhi (e.g. 
Bandiera & Corbelli 2001, for a correct statistical treatment of such measurements). If 
21-cm observations are sensitive enough to Hi column densities where the hydrogen gas is 
mostly ionized in the presence of a UV radiation field, then we could use these data to infer 
g{NH) and the intensity of the ionizing background radiation in the local Universe, which 
is not directly observable. Similarly, the observed sharp drop of the Hi column density in 
the outer disk of two nearby spiral galaxies, M33 and NGC3198, interpreted as an Hl-Hll 
transition zone, made possible an estimate of the UV ionizing flux z — (Corbelli & 
Salpeter 1993; Maloney 1993). The known dark matter content of these galaxies, inferred 
from the gas kinematics, reduced the associated uncertainties. In this paper we shall give 
additional evidence for a UV ionizing background in the local Universe using the shape of 
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the observed f{NHi)- We derive power law fits to giNn) which reproduce f{NHi) quite 
well over many orders of magnitude, after ionization corrections have been applied. We 
underline the observations needed in the future to reduce the ambiguities in modeling the 
ionization corrections and the resulting best fitting values of g{NH). We briefiy discuss two 
different evolutionary scenarios for gaseous structures in the Universe, both compatible with 
the present data. 



2. The Nhi distribution function at z = 

At low redshifts there are two different sets of data which can be used to derive /(A^ij/): 

- Lya absorption lines and Lyman breaks observed in QSOs spectra. The Hubble Space 
Telescope (HST) has provided good quality data for low redshift absorption systems. Penton 
et al. (2000) and ShuU (2001) have analyzed HST/FUSE data for redshifts z < 0.1 and 
Nhi < 10^^ cm~^. Weymann et al. (1998) have analyzed data for N^i < 10^^'^ cm~^ in 
a broader redshift range {z < 1.3) and found a weak evolution in the number of absorbers 
per unit redshift. As Penton et al. (2000) have already pointed out there is good agreement 
between the Nhi distributions derived from the two sets of data and therefore we limit 
ourselves to the Weymann et al. (1998) data, shown in Figure 1(a). For higher column 
densities /{Nhi) at 2; ~ is evaluated from the Lyman limit and Damped Lya data collection 
of Bandiera & Corbelli (2001) over the redshift interval < 2; < 1.3, taking into account the 
mild evolution of Lyman limit systems (CorbeUi et al. 2001). 

- 21-cm maps of nearby galaxies and the Hi mass function. 21-cm observation of the 
Hi gas in nearby galaxies with high resolution and sensitivity provide data on dE{NHi), 
the galaxy differential cross section for a given range of Nhi, averaged over all inclinations. 
This information, acquired for galaxies of different Hi masses and coupled with the Hi mass 
function, ^Mhi), gives fiNni): 

^^L ^{MHi)dT.{NHi,MHi)dMHi 

= (H lA dN 

[Hq/c) dNni 

In order to determine f{NHi) across the region where the hydrogen becomes mostly ionized, 
observations should be sensitive to column densities down to ~ 3 x 10^^ cm~^. Recently the 
Hi Parkes All Sky Survey (HIPASS) together with the Austrafia Telescope Compact Array 
(ATCA) has provided good quality data for (^{Mhi) and for dT,{NHi, Mhi) (Ryan- Weber et 
al. 2001). Fitting the data with a Schechter function (Schechter 1976), the resulting ^{Mhi) 
has a characteristic Hi mass value M* — 1.14 x 10^° M0, a faint end slope a — 1.52, and a 
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normalization 6 — 0.0032 (Kilborn 1999, 2000). Similar values have been obtained recently 
also by Rosenberg & Schneider (2001). We use this mass function in the following analysis 
and consider 6 mass bins of equal width across the range 7.5 < log(M/j//MQ)< 10.5. Ryan- 
Weber et al. (2001) have evaluated dT.{NHi, Mhi) for two Mhi ranges: for lO"^-^ < Mhi < 
10^ M0 and for 10^° < Mhi < 10^°-^ Mq. In the literature there are furthermore two 
galaxies which have been mapped in great detail down to a low surface density: M33 and 
NGC3198. We compute for these galaxies as orientation averages of the face on cross 
sections, obtained from the deconvolution of the observed Nhi radial profiles with the best 
fitted tilted ring model (Begeman 1989; Maloney 1993; CorbeUi & Salpeter 1993; CorbeUi & 
Salucci 2000). The Hi mass of M33 is ~ 2.7 x lO'^ Mq for a distance of 0.84 Mpc and that 
of NGC3198 is ~ 5 X 10^ M0 for a distance of 9.4 Mpc. Data on these two galaxies can 
therefore be used for dTi{NHi, Mhi) in the two mass bins spanning 10^ < Mhi < 10^'^ Mq. 
In order to estimate rfS for 10® < Mhi < 10^ Mq, where we have no data, we have averaged 
the values on the adjacent Mhi bins. /{Nhi), as derived from equation (1), is plotted in 
Figure 1(6). The ATCA data (Ryan- Weber et al. 2001), as weU as the M33 and NGC3198 
data, show the "footprint" of ionization. While in 21-cm emission maps the ionization has 
the effect of giving a sharp drop to the radial decline of Nhi, in the distribution function 
this translates into a plateau where the Hl-Hll transition occurs. When the gas becomes 
mostly ionized, at lower Nhi values, /{Nhi) rises again but it is offset with respect to the 
neutral side and has a slope which depends on how the Hl/Hll ratio varies with Nh (CorbeUi 
et al. 2001). Figure 1(6) shows that this is indeed the behavior of /{Nhi), even though the 
sum over different mass bins may have smoothed somewhat the Hl-Hll transition (if the gas 
vertical stratification depends on the galaxy mass). 



3. From the shape of /{Nhi) to the total gas distribution 

In order to evaluate /{Nhi) over the range 10^^ < Nhi < 6 x 10^° cm~^ in the local 
Universe we select QSOs absorption data for Nhi < 5 x 10^'^ cm~^ and 21-cm emission data 
at higher Nhi- We exclude both the very low column density Lya forest, since it may not be 
associated with denser structures, and the very high column density data, because star and 
molecular gas formation might have reduced the original Hi column density. While simple 
power law fits to /{Nhi) fail we will show in what follows that a power law for the total 
hydrogen distribution can instead well fit all the data displayed in Figure 2. As in CorbeUi 
et al. (2001) we derive a theoretical /{Nhi) after applying the Nh — Nhi conversion relation 
to an assumed power law distribution for the total gas column density, g{NH) — KN^'^. 
The use of 21-cm data eliminates the problem of having large errors in Nhi in the region 
where Lyman limit systems become saturated and therefore we can safely use /{Nhi) binned 
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over small Nhi intervals. We consider two possible Nh — Nhi conversion models for 10^ < 



-(a) the model used by Corbelli et al. (2001) to fit the data on Lyman limit and Damped 
Lya systems, as described in their Section 3, extended to lower column density systems. This 
model is valid for a dark matter confined gas with ionization fractions which increase steadily 
as the total gas column density decreases. 

-(b) the model used by CorbeUi et al. (2001) only for Nh > N^; for Nh < N^ we keep 
constant the H ionization fraction, equal to that reached for Nh — N^. A constant ionization 
fraction for Nh < N^ is achieved if the average gas volume density stays constant as Nh 
decreases below N^. A constant gas volume density might imply either that a constant 
external pressure is confining the gas or that dark matter and turbulence are distributed in 
such a way as to give a gas vertical dispersion which does not depend on the gas column 
density. 

For both (a) and (6) we find the best fitting value of a and of Jl/vo- This is the ratio 
between the intensity of the background flux at 912 A normalized to 10~^^ ergs cm~^ s~^ 
Hz~^ sr~^, and the gas compression factor rjo due to gravity from matter other than the 
gas (e.g. dark matter). The constant rjo adds to the gas self gravity in vertical scale height 
expression (CorbeUi et al. 2001): 



where Mg is the total mass surface density per cm of g diS , Cg IS the thermal sound 
speed and is the value of Cg for T = 10^ K. For model (a) the minimum value of (1-3) 
corresponds to a = 3.32 and to log Jl/tjo = —0.25. The 3-cr uncertainties on a are ±0.4, 
on logJ^/r^Q are ±0.6. An a — 2> power law is therefore compatible with the data. This 
is important since there have been often claims that /(Nhi) for disk hke structures should 
follow an Nff^ power law (Milgrom 1988), but /{Nhi) observed in galaxies strongly deviates 
from this law (Zwaan et al. 1999). This is not surprising since it is giNn) which should 
eventually scale as N^^ for disk like systems, and we have shown that for model (a) the 
present data are compatible with this law. For model (6) the best fit gives a = 1.47 ±0.1, 
log Jl/?7o = —1.85 ± 1, and logA^*/cm~^ = 19.4 ± 0.4. This model fits the data better than 
model (a) with 0.5 being the minimum value. 

The resulting /{Nhi) arc shown in Figure 2 together with the data used for the 
minimization. Both best fitting models reproduce rather well also the high and low column 
density QSO's absorption data, not used for the fit, but shown in the Figure as well. A 



Nhi < 6 X 10- 
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pure power law with no ionization corrections fails since it gives a minimum ^ ~ 10. 

Unfortunately the present data are not sufficient to discriminate between (a) and (6) type 
of model. This is mostly due to the lack of data around Nhi — 10^^ cm~^ and to the 
uncertainties on the conversion of atomic hydrogen gas into other baryonic forms at very 
high column densities. In these two regions (a) and (6) model predictions for f{NHi) differ 
substantially. 

The two best fitting power laws for qINh) imply a quite different degree of ionization 
and therefore different total gas densities in the local Universe. For model (a) most of 
the gas resides in low column density absorbers; H neutral fractions are of order 0.002 for 
Nhi — 10^^ cm~^ and are as low as 10~^ for Nhi — 10^^ cm~^. The total gas density(H+He) 
predicted by this distribution for 10^^ < Nhi < 10^^ gives ilgh^Q{z = 0) ~ 0.01 with a 3a 
range of 0.005 — 0.025. Model (6) instead predicts a lower total density of gas in the same 
range of Nhi: ^ghoiz = 0) ~ lO^^ with a 3a range of 0.001 - 0.005. 

4. Discussion 

We now discuss the best fitting power laws for g^Nu) derived in the previous Sections, 
taking into account the known limits on the UV ionizing background flux at z = 0, the limits 
on the density of baryonic matter, Q^, and the total gas column density distribution obtained 
for QSOs absorbers at 2; ~ 2 — 3 by CorbeUi et al. (2001). Limits on the intensity of the UV 
ionizing background at ^ = rely on several indirect techniques such as those based on the 
detection of a sharp cutoff in the Hi surface density of outer disks of nearby galaxies (CorbeUi 
& Salpeter 1993; Maloney 1993; Dove & Shull 1994), on searches for Ha emission from outer 
disks or from intergalactic Hi clouds (Donahue et al. 1995; Madsen et al. 2001; Weymann 
et al. 2001), on measurements of column densities of metal species in intergalactic Hi clouds 
(Tumlinson et al. 1999). The range Jl — 0.1 — 0.4 is consistent with the above hmits as 
well as with the determination of the UV ionizing background from the QSOs luminosity 
function, even considering possible additional contributions of photons escaping from star 
forming galaxies (Haardt & Madau 1996; Shull et al. 1999). Measurements of the proximity 
effect for 2; < 1 have been reported (Kulkarni & Fall 1993; Pascarelle et al. 2001; Scott et al. 
2001) but the resulting ionizing background intensity, Jl — 0.01 —4, is highly uncertain due 
to possible clustering of the Ly« forest around QSOs and to a possible evolution over the used 
z range. The total density of baryonic matter can be inferred from deuterium abundances, 
fife/i^o ~ 0.02 (Buries et al. 2001; O'Meara et al. 2001). For a Universe with = 0.7 and 
JIm = 0.3 this Qft value limits the power law index of the total gas column density distribution 
derived by CorbeUi et al. (2001) at higher redshifts to values a < 3. The extrapolation of 
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the same power law towards the Lya forest requires however that ionization conditions of 
the gas stay constant for Nhi < 10^^ cm~^ at these redshifts i.e. for Nh ^ 2.5 x 10^^ cm~^. 
In the foUowing we shah assume that the gas distribution at z ~ 2 — 3 is well described by 
the Best-Fit 1 model of CorbeUi et al. (2001) {a = 2.7, r]o = 12.5, log J^/t^q = -0.35, see 
their Table 1). 

If today the gas is distributed according to the best fitting model (a), with a ~ 3.3, the 
limits on and on at 2; = force the gas compression factor to be r^o < 3. According 
to this model the following evolutionary scenario might hold: because changes in time of 
the ionization conditions of the gas arc mild, the average gas compression factor should 
decrease with time in order to balance the decrease of Jl- If the dark matter is confining 
and compressing the gas then today's systems have on average a lower dark matter density 
than systems at earlier times. The alternative would be to invoke a steadily increasing 
turbulence which lowers the effective compression 770, or additional ionization sources in 
the surroundings of gaseous clouds. The power law index for the total gas column density 
distribution from z ~ 2.5 to z ~ becomes slightly steeper and now it is compatible with 
a disk-like geometry. For model (a) the gas associated with the Lya forest clouds is, even 
today, the dominant reservoir of baryonic matter in the Universe. 

The best fitting model (6) aX z ^ implies a constant hydrogen neutral fraction, of 
order 0.05, for Nhi ^ 10^^ cm^^. The power law index for g{NH) in this case is a ^ 1.5 
and Jl/i1o — 0.01. In order to explain the flattening of the power law index for giNn) 
with time, we invoke a scenario in which low column density structures collapse towards 
the center or merge to form higher column density structures. If rjo is constant through the 
evolution of the Universe, the required ionizing background radiation at 2; = would equal 
that estimated by Haardt & Madau (1996) for QSOs as primary sources of ionizing photons. 
If the external pressure is responsible for confining the gas and keeping constant the ionizing 
fraction at low cohimn densities, the inferred value of Pg^t is also constant through redshifts 
and is Pe^t /{cm-'^ K) ~ r/o ~ 10 (Charlton et al. 1994; CorbeUi & Salpeter 1995). This 
excludes that Pext is dominated by a diffuse and pervasive intergalactic medium since this is 
expected to decrease with time. However, the merging or the collapse of low column density 
structures is likely to trigger star formation and outflows of hot gas into the intergalactic 
medium; this hot gas at temperatures T~ 10^-10^ K would be the hkely confining medium 
of low column density Lya clouds and it would contribute mostly to the density of baryonic 
matter in the Local Universe (Fukugita et al. 1998; Cen & Ostriker 1999; Dave et al. 2001; 
Tripp 2001). 

In this paper we have shown that ionizing radiation affects the neutral gas distribution 
in the local Universe and that any attempt to fit /{Nhi) with a simple power law across 
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the Hl-Hll transition region (10^^ < Nhi < 10^° cm~^) can lead to misleading results. A 
power law distribution for the total gas column density instead fits the data reasonably well, 
from the Lya forest to the high gas column density in galaxies. Additional 21-cm data, as 
well as new Lyman limit searches at low redshifts, are needed to implement the knowledge 
of /{Nhi) in this delicate region and draw a more definite conclusion on the slope of the 
total gas column density distribution. Furthermore, if 21-cm data are analyzed in detail 
to constrain the gas and dark matter distribution in outer disks, one can assign a more 
definite value to the gas compression and picture a unique evolutionary scenario for the gas 
distribution and the UV ionizing background in the Universe. 

We are grateful to the anonymous referee for a careful reading of the original manuscript. 
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Fig. 1. — (a) f{NHi) from QSO Lya and Lyman continuum absorption data at low redshifts. 
Lya forest data from Weymann et al. (1998) are shown as filled triangles while data for 
Lyman limit and Damped Lya systems from our collection are shown as open dots, (b) 
/{Nhi) from 21-cm data at 2; = (see text for details). 
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Fig. 2.— The two theoretical J^Nhi) which best fit the data for 10^^ < Nhi < 6 x 10^° cm-^, 
shown in the Figure with simple crosses. The continuous line is derived from a total gas 
distribution function with a = 3.32 using the ionization model (a) with log J/r/o = —0.25. 
The dashed line is derived from a total gas distribution function with a — 1.47 using the 
ionization model (b) with log J/770 = —1-9 and — 2.5 x 10^^ cm~^. The high and low 
column density data shown with filled dots have not been used for the fit. 



